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(54) Title: METHOD AND DEVICE FOR IDENTIFYING A BIOLOGICAL SAMPLE 

(57) Abstract: Methods, apparatus and systems for identifying a biological sample that generate a data set indicative of the com- 
position of the biological sample are provided. In a particular example, the data set is DNA spectrometry data received from a mass 
spectrometer. The data set is denoised, and a baseline is deleted. Since possible compositions of the biological sample may be 
known, expected peak areas may be determined. Using the expected peak areas, a residual baseline is generated to further correct the 
data set. Probable peaks are then identifiable in the corrected data set, which are used to identify the composition of the biological 
sample. In a disclosed example, statistical methods are employed to deterrnine the probability that a probable peak is an actual peak, 
not an actual peak, or that the data are too inconclusive to call. 
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METHOD AND DEVICE FOR IDENTIFYING A BIOLOGICAL SAMPLE 
RELATED APPLICATIONS 

Benefit of priority is claimed to U.S. application Serial No. 09/663,968, 
filed September 19, 2000, entitled "METHOD AND DEVICE FOR IDENTIFYING A 
5 BIOLOGICAL SAMPLE" to Ping Yip. Where permitted the subject matter of this 
application is incorporated in its entirety by reference. 

This application is related to U.S. application Serial No. 09/285,481, filed 
April 2, 1999, entitled "AUTOMATED PROCESS LINE", to Hubert Koster, Ping 
Yip, Jhobe Steadman, Dirk Reuter and Richard MacDonald. Where permitted the 
10 subject matter of this application is incorporated in its entirety by reference. 
HELD OF THE INVENTION 

The present invention is in the field of biological identification. More 
specifically, processes and systems for identifying a biological sample by 
analyzing information received from a test instrument are provided. 
1 5 BACKGROUND OF THE INVENTION 

Advances in the field of genomics are leading to the discovery of new 
and valuable information regarding genetic processes and relationships. This 
newly illuminated genetic information is revolutionizing the way medical 
therapies are advanced, tested, and delivered. As more information is gathered, 
20 genetic analysis has the potential to play an integral and central role in 

developing and delivering medical advancements that will significantly enhance 
the quality of life. 

With the increasing importance and reliance on genetic information, the 
accurate and reliable collection and processing of genetic data is critical. 

25 However, conventional known systems for collecting and processing genetic or 
DNA data are inadequate to support the informational needs of the genomics 
community. For example, known DNA collection systems often require 
substantial human intervention, which undesirably risks inaccuracies associated 
with human intervention. Further, the slow pace of such a manual task severely 

30 limits the quantity of data that can be collected in a given period of time, which 
slows needed medical advancements and adds substantially to the cost of data 
collection. 
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In a particularly exciting area of genomics, the identification and 
classification of minute variations in human DNA has been linked with 
fundamental treatment or medical advice for a specific individual. For example, 
the variations are a strong indication of predisposition for a particular disease, 
5 drug tolerance, and drug efficiency. The most promising of these minute 

variations are commonly referred to as Single Nucleotide Polymorphisms (SNPs), 
which relate to a single base-pair change between a first subject and a second 
subject. By accurately and fully identifying such SNPs, a health care provider 
would have a powerful indication of a person's likelihood of succumbing to a 
10 particular disease, which drugs will be most effective for that person, and what 
drug treatment plan will be most beneficial. Armed with such knowledge, the 
health care provider can assist a person in lowering other risk factors for high- 
susceptibility diseases. Further, the health care provider can confidently select 
appropriate drug therapies, a process which is now an iterative, hit or miss 
1 5 process where different drugs and treatment schedules are tried until an 

effective one is found. Not only is this a waste of limited medical resources, but 
the time lost in finding an effective therapy can have serious medical 
consequences for the patient. 

In order to fully benefit from the use of SNP data, vast quantities of DNA 
20 data must be collected, compared, and analyzed. For example, collecting and 
identifying the SNP profile for a single human subject requires the collection, 
identification, and classification of thousands, even tens of thousands of DNA 
samples. Further, the analysis of the resulting DNA data must be carried out 
with precision. In making a genetic call, where a composition of a biological 
25 sample is identified, any error in the call may result in detrimentally affecting the 
medical advice or treatment given to a patient. 

Conventional, known systems and processes for collecting and analyzing 
DNA data are Inadequate to timely and efficiently implement a widespread 
medical program benefiting from SNP information. For example, many known 
30 DNA analysis techniques require the use of an operator or technician to monitor 
and review the DNA data. An operator, even with sufficient training and 
substantial experience, is still likely to occasionally make a classification error. 
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For example, the operator may incorrectly identify a base-pair, leading to that 
patient receiving faulty SNP profile. Alternatively, the operator may view the 
data and decide that the data do not clearly Identify any particular base pair. 
Although such a "no call" may be warranted, it is likely that the operator will 
5 make "no-cal!" decisions when the data actually support a valid call. In such a 
manner, the opportunity to more fully profile the patient is lost. 

Thus, there exists a need for systems, apparatus and processes to 
efficiently and accurately collect and analyze data, such as DNA data. 
Thererfore it is an object herein to provide such systems., apparatus and 
10 processes, it Is an object herein to provide an apparatus and process for 
accurately identifying genetic information. It Is another object to provide 
processes and apparatus for extracting genetic information from genetic data in 
a highly automated manner. To overcome the deficiencies In the known 
conventional systems, a method and apparatus for Identifying a biological 
1 5 sample are provided. 

SUMMARY OF THE INVENTION 

A method and system for Identifying a biological sample that generates a 
data set indicative of the composition of the biological sample are provided. In a 
particular example, the data set is DNA spectrometry data received from a mass 
20 spectrometer. The data set Is denoised, and a baseline is deleted. Since 

possible compositions of the biological sample may be known, expected peak 
areas may be determined. Using the expected peak areas, a residual baseline is 
generated to further correct the data set. Probable peaks are then Identifiable fn 
the corrected data set, Which are used to identify the composition of the 
25 biological sample. In a disclosed example, statistical methods are employed to 
determine the probability that a probable peak is an actual peak, not an actual 
peak, or that the data are too Inconclusive to call. 

Advantageously, the method and system for identifying a biological 
samplb accurately makes composition calls In a highly automated manner. In 
30 such a manner, complete SNP profile information, for example, may be collected 
efficiently. More importantly, the collected data are analyzed with highly 
accurate results. For example, when a particular composition is called, the resuit 
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may be relied upon with great confidence. Such confidence is provided by the 
robust computational process employed and the highly automatic method of 
collecting, processing, and analyzing the data set. 

These and other features and advantages of the present invention will be 
5 appreciated from review of the following detailed description of the invention, 
along with the accompanying figures in which like reference numerals refer to 
like parts throughout. 
BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram showing a system provided herein; 
10 FIG. 2 is a flowchart of a method of identifying a biological sample 

provided herein; 

FIG. 3 is a graphical representation of data from a mass spectrometer; 
FIG. 4 is a diagram of wavelet transformation of mass spectrometry data; 
FIG. 5 is a graphical representation of wavelet stage 0 hi data; 
15 FIG. 6 is a graphical representation of stage 0 noise profile; 

FIG. 7 is a graphical representation of generating stage noise standard 
deviations; 

FIG. 8 is a graphical representation of applying a threshold to data stages; 

FIG. 9 is a graphical representation of a sparse data set; 
20 FIG. 10 is a formula for signal shifting; 

FIG. 1 1 is a graphical representation of a wavelet transformation of a 
denoised and shifted signal; 

FIG. 12 is a graphical representation of a denoised and shifted signal; 

FIG. 13 is a graphical representation of removing peak sections; 
25 FIG. 14 is a graphical representation of generating a peak free signal; 

FIG. 1 5 is a block diagram of a method of generating a baseline 
correction; 

FIG. 16 is a graphical representation of a baseline and signal; 
FIG. 17 is a graphical representation of a signal with baseline removed; 
30 FIG. 1 8 is a table showing compressed data; 

FIG. 1 9 is a flowchart of method for compressing data; 
FIG. 20 is a graphical representation of mass shifting; 
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FIG. 21 is a graphical representation of determining peak width; 
RG. 22 is a graphical representation of removing peaks; 
FIG. 23 is a graphical representation of a signal with peeks removed; 
FIG. 24 is a graphical representation of a residual baseline; 
5 FIG. 25 is a graphical representation of a signal with residual baseline 

removed; 

FIG. 26 is a graphical representation of determining peak height; 
FIG. 27 is a graphical representation of determining signal-to-nolse for 
each peak; 

10 FIG. 28 is a graphical representation of determining a residual error for 

each peak; 

RG. 29 Is a graphical representation of peak probabilities; 
RG. 30 is a graphical representation of applying an allelic ratio to peak 
probability; 

15 FIG. 31 Is a graphical representation of determining peak probability; 

FIG. 32 is a graphical representation of calling a genotype; and 
FIG. 33 is a flowchart showing a statistical procedure for calling a 
genotype. 

DETAILED DESCRIPTION OF THE INVENTION 

20 Provided herein are a method and device for identifying a biological 

sample. Referring now to FIG, 1, an apparatus 10 for identifying a biological 
sample is disclosed. The apparatus 10 for identifying a biological sample 
generally comprises a mass spectrometer 1 6 communicating with a computing 
device 20. In a preferred embodiment, the mass spectrometer may be a MALDI- 

25 TOF <T>ass spectrometer manufactured by Bruker-Franzen Analytik GmbH; 
however, it will be appreciated that other mass spectrometers can be 
substituted. The computing device 20 is preferably a general purpose computing 
device. However, it will be appreciated that the computing device could be 
alternatively configured; for example, it may be integrated with the mass 

30 spectrometer or could be part of a computer in a larger network system. 

The apparatus 10 for identifying a biological sample may operate as an 
automated identification system having a robot 25 with a robotic arm 27 
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configured to deliver a sample plate 29 into a receiving area 31 of the mass 
spectrometer 15. In such a manner, the sample to be identified may be placed 
on the plate 29 and automatically received into the mass spectrometer 15. The 
biological sample is then processed in the mass spectrometer to generate data 
5 indicative of the mass of DNA fragments into biological sample. These data may 
be sent directly to computing device 20, or may have some preprocessing or 
filtering performed within the mass spectrometer. In a preferred embodiment, 
the mass spectrometer 1 5 transmits unprocessed and unf iltered mass 
spectrometry data to the computing device 20. However, it will be appreciated 
10 that the analysis in the computing device may be adjusted to accommodate 
preprocessing or filtering performed within the mass spectrometer. 

Referring now to FIG. 2, a general method 35 for identifying a biological 
sample is shown. In method 35, data are received into a computing device from 
a test instrument in block 40. Preferably the data are received in a raw, 
15 unprocessed and unfiltered form, but alternatively may have some form of 

filtering or processing applied. The test instrument of a preferred embodiment is 
a mass spectrometer as described above. However, it«will be appreciated that 
other test instruments could be substituted for the mass spectrometer. 

The data generated by the test instrument, and in particular the mass 
20 spectrometer, include information indicative of the identification of the biological 
sample. More specifically, the data are indicative of the DNA composition of the 
biological sample. Typically, mass spectrometry data gathered from DNA 
samples obtained from DNA amplification techniques are noisier than, for 
example, those from typical protein samples. This is due in part because protein 
25 samples are more readily prepared in more abundance, and protein samples are 
more easily ionizable as compared to DNA samples. Accordingly, conventional 
mass spectrometer data analysis techniques are generally ineffective for DNA 
analysis of a biological sample. 

To improve the analysis capability so that DNA composition data can be 
30 more readily discerned, a preferred embodiment uses wavelet technology for 
analyzing the DNA mass spectrometry data. Wavelets are an analytical tool for 
signal processing, numerical analysis, and mathematical modeling. Wavelet 



WO 02/25567 



PCTYUS01/29290 



\ 

-7- 

technology provides a basic expansion function which is applied to a data set. 
Using wavelet decomposition, the data set can be simultaneously analyzed in 
both the time and frequency domains. Wavelet transformation is the technique 
of choice in the analysis of data that exhibit complicated time (mass) and 
5 frequency domain information, such as MALDI-TOF DNA data. Wavelet 

transforms as described herein have superior denoising properties as compared 
to conventional Fourier analysis techniques. Wavelet transformation has proven 
to be particularly effective in interpreting the inherently noisy MALDI-TOF 
spectra of DNA samples. In using wavelets, a "small wave" or "scaling 
10 function" is used to transform a data set into stages, with each stage 
representing a frequency component in the data set. Using wavelet 
transformation, mass spectrometry data can be processed, filtered, and analyzed 
with sufficient discrimination to be useful for identification of the DNA 
composition for a biological sample. 
15 Referring again to FIG. 2, the data received in block 40 are denoised in 

block 45. The denoised data then has a baseline correction applied in block 50. 
A baseline correction is generally necessary as data coming from the test 
instrument, in particular a mass spectrometer instrument, has data arranged in a 
generally exponentially decaying manner. This generally exponential decaying 
20 arrangement is not due to the composition of the biological sample, but is a 
result of the physical properties and characteristics of the test instrument and 
other chemicals involved in DNA sample preparation. Accordingly, baseline 
correction substantially corrects the data to remove a component of the data 
attributable to the test system and sample preparation characteristics. 
25 After denoising in block 45 and the baseline correction in block 50, a 

signal remains which is generally indicative of the composition of the biological 
sample. However, due to the extraordinary discrimination required for analyzing 
the DNA composition of the biological sample, the composition is not readily 
apparent from the denoised and corrected signal. For example, although the 
30 signal may include peak areas, it is not yet clear whether these "putative" peaks 
actually represent a DNA composition, or whether the putative peaks are result 
of a systemic or chemical aberration. Further, any call of the composition of the 
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biological sample would have a probability of error which would be unacceptable 
for clinical or therapeutic purposes. In such critical situations, there needs to be 
a high degree of certainty that any call or identification of the sample Is 
aacurate. Therefore, additional data processing and interpretation are necessary 
5 before the sample can be accurately and confidently identified. 

Since the quantity of data resulting from each mass spectrometry test Is 
typically thousands of data points, and an automated system may be set to 
perform hundreds or even thousands of tests per hour, the quantity of mass 
spectrometry data generated is enormous. To facilitate efficient transmission 

10 and storage of the mass spectrometry data, block 85 shows that the denolsed 
and baseline corrected data are compressed. 

In a preferred embodiment, the biological sample is selected and 
processed to have only a limited range of possible compositions. Accordingly, It 
is therefore known where peaks indicating composition should be located, if 

16 present. Taking advantage of knowing the location of these expected peaks, in 
block 60 the method 35 matches putative peaks in the processed signal to the 
location of the expected peaks. In such a manner, the probabilitY of each 
putative peak in the data being an actual peak Indicative of the composition of 
the biological sample can be determined. Once the probability of each peak is 

20 determined In block 60, then in block 65 the method 36 statistically determines 
the composition of the biological sample and determines if confidence is high 
enough to calling a genotype. 

Referring again to block 40, data are received from the test instrument, 
which is preferably a mass spectrometer. In a specific illustration, FIG. 3 shows 

25 an example of data from a mass spectrometer. The mess spectrometer data 70 
generally comprises data points distributed along an x-axis and a y-axis. The x- 
axis represents the mass of particles detected, while the y-axis represents a 
numerical concentration of the particles. As can be seen in FIG. 3, the mass 
spectrometry data 70 is generally exponentially decaying with data at the left 

30 end of the x-axis generally decaying in an exponential manner toward data at the 
heavier end of the x-axis. However, the general exponential presentation of the 
data is not Indicative of the composition of the biological sample, but is more 
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reflective of systematic error and characteristics. Further, as described above 
and illustrated in FIG. 3, considerable noise exists in the mass spectrometry DNA 
data 70. 

Referring again to block 45, where the raw data received in block 40 is 
5 denoised, the denoising process will be described in more detail. As illustrated 
in FIG. 2, the denoising process generally entails 1) performing a wavelet 
transformation on the raw data to decompose the raw data into wavelet stage 
coefficients; 2) generating a noise profile from the highest stage of wavelet 
coefficients; and 3) applying a scaled noise profile to other stages in the wavelet 
10 transformation. Each step of the denoising process is further described below. 

Referring now to FIG. 4, the wavelet transformation of the raw mass 
spectrometry data is generally diagramed. Using wavelet transformation 
techniques, the mass spectrometry data 70 is sequentially transformed into 
stages. In each stage the data is represented in a high stage and a low stage, 
15 with the low stage acting as the input to the next sequential stage. For 

example, the mass spectrometry data 70 is transformed into stage 0 high data 
82 and stage 0 low data 83. The stage 0 low data 83 is then used as an input 
to the next level transformation to generate stage 1 high data 84 and stage 1 
low data 85. In a similar manner, the stage 1 low data 85 is used as an input to 
20 be transformed into stage 2 high data 86 and stage 2 low data 87. The 

transformation is continued until no more useful information can be derived by 
further wavelet transformation . For example, in the preferred embodiment a 24- 
point wavelet is used. More particularly a wavelet commonly referred to as the 
Daubechies 24 is used to decompose the raw data. However, it will be 
25 appreciated that other wavelets can be used for the wavelet transformation. 

Since each stage in a wavelet transformation has one-half the data points of the 
previous stage, the wavelet transformation can be continued until the stage n 
low data 89 has around 50 points. Accordingly, the stage n high 88 would 
contain about 100 data points. Since the preferred wavelet is 24 points long, 
30 little data or information can be derived by continuing the wavelet transformation 
on a data set of around 50 points. 
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FK3. 5 shows an example of stage 0 high data 95. Since stage 0 high 
data 95 is generally indicative of the highest frequencies in the mass 
spectrometry data, stage 0 high data 95 wil) closely relate to the quairtity of 
high frequency noise (n the mass spectrometry data. In FIG, 6, an exponential 
5 fitting formula has been applied to the stage 0 high data 96 to generate a stage 
o noise profile 97. In particular, the exponential fitting formula is In the format 
Aq + A, EXP {-A 2 mJ. It will he appreciated that other exponential fitting 
formulas or other types of curve fits may be used. 

Referring now to FIG. 7 # noise profiles for the other high stages are 

10 determined. Since the later data points In each stage will likely be representative 
of the level of noise In each stage, only the later data points in each stage are 
used to generate a standard deviation figure that is representative of the noise 
content in that particular stage. More particularly, in generating the noise profile 
for each remaining stage, only the last five percent of the data points in each 

15 stage are analyzed to determine a standard deviation number. It will be 

appreciated that other numbers of points or alternative methods could be used to 
generate such a standard deviation figure. 

The standard deviation number for each stage is used with the stage 0 
noise profile (the exponential curve) 97 to generate a scaled noise profile for 

20 each stage. For example, FIG. 7 shows that stage 1 high data 98 has stage 1 
high data 1 03 with the last five percent of the data points represented by area 
99. The points in area 99 are evaluated to determine a standard deviation 
number indicative of the noise content in stage 1 high data 103. The standard 
deviation number is then used with the stage 0 noise profile 97 to generate a 

25 stage 1 noise profile. 

In a simitar manner, stage 2 high 100 has stage 2 high data 104 with the 
last five percent of points represented by area 101 , The data points in area 101 
are then used to calculate a standard deviation number which is then used to 
scale the stage 0 noise profile 97 to generate a noise profile for stage 2 data. 

30 This same process is continued for each of the stage high data as shown by the 
stage n high 105. For stage n high 105, stage n high data 108 has the last five 
percent of data points indicated in area 106. The data points in area 106 are 
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used to determine a standard deviation number for stage n. The stage n 
standard deviation number is then used with the stage 0 noise profile 97 to 
generate a noise profile for stage n. Accordingly, each of the high data stages 
has a noise profile. 

5 FIG. 8 shows how the noise profile is applied to the data in each stage. 

Generally, the noise profile is used to generate a threshold which is applied to 
the data in each stage. Since the noise profile is already scaled to adjust for the 
noise content of each stage, calculating a threshold permits further adjustment 
to tune the quantity of noise removed. Wavelet coefficients below the threshold 
10 are ignored while those above the threshold are retained. Accordingly, the 
remaining data has a substantial portion of the noise content removed. 

Due to the characteristics of wavelet transformation, the lower stages, 
such as stage 0 and 1 , will have more noise content than the later stages such 
as stage 2 or stage n. Indeed, stage n low data is likely to have little noise at 

15 all. Therefore, in a preferred embodiment the noise profiles are applied more 
aggressively in the lower stages and less aggressively in the later stages. For 
example, FIG. 8 shows that stage 0 high threshold is determined by multiplying 
the stage 0 noise profile by a factor of four. In such a manner, significant 
numbers of data points in stage 0 high data 95 will be below the threshold and 

20 therefore eliminated. Stage 1 high threshold 1 12 is set at two times the noise 
profile for the stage 1 high data, and stage 2 high threshold 1 14 is set equal to 
the noise profile for stage 2 high. Following this geometric progression, stage n 
high threshold 1 1 6 is therefore determined by scaling the noise profile for each 
respective stage n high by a factor equal to (1/2*- 2 ). It will be appreciated that 

2B other factors may be applied to scale the noise profile for each stage. For 
example, the noise profile may be scaled more or less aggressively to 
accommodate specific systemic characteristics or sample compositions.. As 
indicated above, stage n low data does not have a noise profile applied as stage 
n low data 1 1 8 is assumed to have little or no noise content. After the scaled 

30 noise profiles have been applied to each high data stage, the mass spectrometry 
data 70 has been denoised and is ready for further processing. A wavelet 
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transf ormation of the denoised signal results in the sparse data set 1 20 as 
shown in FIG. 9- 

Referring again to FIG. 2, the mass spectrometry data received in block 
40 has been denoised in block 45 and is now passed to block 50 for baseline 
5 correction. Before performing baseline correction, the artifacts introduced by the 
wavelet transformation procedure are preferably removed. Wavelet 
transformation results vary slightly depending upon which point of the wavelet is 
used as a starting point. For example, the preferred embodiment uses the 24- 
point Daubechies-24 wavelet. By starting the transformation at the 0 point of 
10 the wavelet, a slightly different result will be obtained than if starting at points 1 
or 2 of the wavelet. Therefore, the denoised data is transformed using every 
available possible starting point, with the results averaged to determine a final 
denoised and shifted signal. For example, FIG. 10 shows that the wavelet 
coefficient is applied 24 different times and then the results averaged to 
15 generate the final data set. It will be appreciated that other techniques may be 
used to accommodate the slight error introduced due to wavelet shifting. 

The formula is generally indicated in FIG. 10. Once the signal has been 
denoised and shifted, a denoised and shifted signal 1 30 is generated as shown 
in FIG. 12. FIG. 11 shows an example of the wavelet coefficient 135 data set 
20 from the denoised and shifted signal 130. 

FIG. 13 shows that putative peak areas 145, 147, and 149 are located in 
the denoised and shifted signal 150. The putative peak areas are systematically 
identified by taking a moving average along the signal 150 and identifying 
sections of the signal 1 50 which exceed a threshold related to the moving 
25 average. It will be appreciated that other methods can be used to identify 
putative peak areas in the signal 1 50. 

Putative peak areas 145, 147 and 149 are removed from the signal 150. 
to create a peak-free signal 1 55 as shown in FIG. 1 4. The peak-free signal 1 55 
is further analyzed to identify remaining minimum values 1 57, and the remaining 
30 minimum values 157 are connected to generate the peak-free signal 155. 

FIG. 15 shows a process of using the peak-free signal 155 to generate a 
baseline 170 as shown in FIG. 16. As shown in block 162, a wavelet 
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transformation is performed on the peak-free signal 155. All the stages from the 
wavelet transformation are eliminated in block 1 64 except for the n low stage. 
The n low stage will generally indicate the lowest frequency component of the 
peak-free signal 155 and therefore will generally indicate the system exponential 
5 characteristics. Block 166 shows that a signal is reconstructed from the n low 
coefficients and the baseline signal 170 is generated in block 168. 

FIG. 16 shows a denoised and shifted data signal 172 positioned adjacent 
a correction baseline 170. The baseline correction 170 is subtracted from the 
denoised and shifted signal 172 to generate a signal 175 having a baseline 
10 correction applied as shown in FIG 17. Although such a denoised, shifted, and 
corrected signal is sufficient for most identification purposes, the putative peaks 
in signal 175 are not identifiable with sufficient accuracy or confidence to call 
the DNA composition of a biological sample. 

Referring again to FIG. 2, the data from the baseline correction 50 is now 
15 compressed in block 55; the compression technique used in a preferred 

embodiment is detailed in FIG. 18. In FIG. 18 the data in the baseline corrected 
data are presented in an array format 182 with x-axis points 183 having an 
associated data value 184. The x-axis is indexed by the non-zero wavelet 
coefficients, and the associated value is the value of the wavelet coefficient. In 
20 the illustrated data example in table 182, the maximum value 184 is indicated to 
be 1000. Although a particularly advantageous compression technique for mass 
spectrometry data is shown, it will be appreciated that other compression 
techniques can be used. Although not preferred, the data may also be stored 
without compression. 
25 In compressing the data according to a preferred embodiment, an 

intermediate format 186 is generated. The intermediate format 186 generally 
comprises a real number having a whole number portion 188 and a decimal 
portion 190. The whole number portion is the x-axis point 1 83 while the 
decimal portion is the value data 184 divided by the maximum data value. For 
30 example, in the data 182 a data value *25 n is indicated at x-axis point n 100 n . 
The intermediate value for this data point would be "100.025°. 
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From the intermediate compressed data 1 86 the final compressed data 
195 is generated. The first point of the intermediate data file becomes the 
starting point for the compressed data. Thereafter each data point in the 
compressed data 195 is calculated as follows: the whole number portion (left of 
5 the decimal) is replaced by the difference between the current and the last whole 
number. The remainder (right of the decimal) remains intact. For example, the 
starting point of the compressed data 1 95 is shown to be the same as the 
intermediate data point which is " 100.025". The comparison between the first 
intermediate data point "100.025" and the second intermediate data point 

10 "150,220" is "50.220". Therefore, "50.220" becomes the second point of the 
compressed data 195. In a similar manner, the second intermediate point is 
"150.220" and the third intermediate data point is "500.0001". Therefore, the 
third compressed data becomes "350.000". The calculation for determining 
compressed data points is continued until the entire array of data points is 

15 converted to a single array of real numbers. 

FIG. 1 9 generally describes the method of compressing mass 
spectrometry data, showing that the data file in block 201 is presented as an 
array of coefficients in block 202. The data starting point and maximum is 
determined as shown in block 203, and the intermediate real numbers are 

20 calculated in block 204 as described above. With the intermediate data points 
generated, the compressed data is generated in block 205. The described 
compression method is highly advantageous and efficient for compressing data 
sets such as a processed data set from a mass spectrometry instrument. The 
method is particularly useful for data, such as mass spectrometry data, that use 

25 large numbers and have been processed to have occasional lengthy gaps in x- 
axis data. Accordingly, an x-y data array for processed mass spectrometry data 
may be stored with an effective compression rate of 10x or more. Although the 
compression technique is applied to mass spectrometry data, it will be 
appreciated with the method may also advantageously be applied to other data 

30 sets. 

Referring again to FIG. 2, peak heights are now determined in block 60. 
The first step in determining peak height is illustrated in FIG. 20 where the signal 
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210 is shifted left or right to correspond with the position of expected peaks. 
As the set of possible compositions in the biological sample is known before the 
mass spectrometry data is generated, the possible positioning of expected peaks 
is already known. These possible peaks are referred to as expected peaks, such 
5 as expected peaks 212, 214, and 216. Due to calibration or other errors in the 
test instrument data, the entire signal may be shifted left or right from its actual 
position; therefore, putative peaks located in Hie signal, such as putative peaks 
218, 222, and 224 may be compared to the expected peaks 212, 214, and 216, 
respectively. The entire signal is then shifted such that the putative peaks align 

10 more closely with the expected peaks. 

Once the putative peaks have been shifted to match expected peaks, the 
strongest putative peak is identified in FIG. 21. In a preferred embodiment, the 
strongest peak is calculated as a combination of analyzing both the overall peak 
height and area beneath the peak. For example, a moderately high but wide 

15 peak would be stronger than a very high peak that is extremely narrow. With 
the strongest putative peak identified, such as putative peak 225, a Gaussian 
228 curve is fit to the peak 225. Once the Gaussian is fit, the width (W) of the 
Gaussian is determined and will be used as the peak width for future 
calculations. 

20 As generally addressed above, the denoised, shifted, and baseline- 

corrected signal is not sufficiently processed for confidently calling the DNA 
composition of the biological sample. For example, although the baseline has 
< generally been removed, there are still residual baseline effects present. These 
residual baseline effects are therefore removed to increase the accuracy and 

25 confidence in making identifications. 

To remove the residual baseline effects, FIG. 22 shows that the putative 
peaks 218, 222, and 224 are removed from the baseline corrected. signal. The 
peaks are removed by identifying a center line 230, 232, and 234 of the 
putative peaks 218, 222, and 224, respectively and removing an area both to 

30 the left and to the right of the identified center line. For each putative peak, an 
area equal to twice the width (W) of the Gaussian is removed from the left of 
the center line, while an area equivalent to 50 daltons is removed from the right 
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of the center line. It has been found that the area representing 50 daltons is 
adequate to sufficiently remove the effect of salt adducts which may be 
associated with an actual peak. Such adducts appear to the right of an actual 
peak and are a natural effect from the chemistry involved in acquiring a mass 
5 spectrum. Although a 50 Oalton buffer has been selected, it will be appreciated 
that other ranges or methods can be used to reduce or eliminate adduct effects. 

The peaks are removed and remaining minima 247 located as shown in 
FIG. 23 with the minima 247 connected to create signal 245. A quartic 
polynomial is applied to signal 245 to generate a residual baseline 250 as shown 
10 in FIG. 24. The residual baseline 250 is subtracted from the signal 225 to 
generate the final signal 255 as indicated in FIG. 25. Although the residual 
baseline is the result of a quartic fit to signal 245, it will be appreciated that 
other techniques can be used to smooth or fit the residual baseline. 

To determine peak height as shown in FIG. 26, a Gaussian such as 
15 Gaussian 266, 268, arid 270 is fit to each of the peaks, such as peaks 260, 
262, and 264, respectively. Accordingly, the height of the Gaussian is 
determined as height 272, 274, and 276. Once the height of each Gaussian 
peak is determined, then the method of identifying a biological compound 35 can 
move into the genotyping phase 65 as shown in FIG. 2. 
20 An indication of the confidence that each putative peak is an actual peak 

can be discerned by calculating a signal-to-noise ratio for each putative peak. 
Accordingly, putative peaks with a strong signal-to-noise ratio are generally more 
likely to be an actual peak than a putative peak with a lower signal-to-noise 
ratio. As described above and shown in FIG. 27, the height of each peak, such 
25 as height 272, 274, and 276, is determined for each peak, with the height being 
an indicator of signal strength for each peak. The noise profile, such as noise 
profile 97, is extrapolated into noise profile 280 across the identified peaks. At 
the center line of each of the peaks, a noise value is determined, such as noise 
value 282, 283, and 284. With a signal value and a noise value generated, 
30 signal-to-noise ratios can be calculated for each peak. For example, the signal- 
to-noise ratio for the first peak in FIG. 27 would be calculated as signal value 
272 divided by noise value 282, and in a similar manner the signal-to-noise ratio 
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of the middle peak in FIG. 27 would be determined as signal 274 divided by 
noise value 283. 

Although the signal-to-noise ratio is generally a useful indicator of the 
presence of an actual peak, further processing has been found to increase the 
5 confidence by which a sample can be identified. For example, the signal-to- 
noise ratio for each peak in the preferred embodiment is preferably adjusted by 
the goodness of fit between a Gaussian and each putative peak. It is a 
characteristic of a mass spectrometer that sample material is detected in a 
manner that generally complies with a normal distribution. Accordingly, greater 
10 confidence will be associated with a putative signal having a Gaussian shape 
than a signal that has a less normal distribution. The error resulting from having 
a non-Gaussian shape can be referred to as a "residual error". 

Referring to FIG. 28, a residual error is calculated by taking a root mean 
square calculation between the Gaussian 293 and the putative peak 290 in the 
15 data signal. The calculation is performed on data within one width on either side 
of a center line of the Gaussian. The residual error is calculated as: 

20 N 

where G is the Gaussian signal value, R is the putative 
peak value, and N is the number of points from -W to + W. The calculated 
25 residual error is used to generate an adjusted signal-to-noise ratio, as described 
below. 

An adjusted signal noise ratio is calculated for each putative peak using 
the formula (S/N) * EXP K1-W , where S/N is the signal-to-noise ratio, and R is the 
residual error determined above. Although the preferred embodiment calculates 
30 an adjusted signal-to-noise ratio using a residual error for each peak, it will be 

appreciated that other techniques can be used to account for the goodness of fit 
between the Gaussian and the actual signal. 

Referring now to FIG. 29, a probability is determined that a putative peak 
is an actual peak. In making the determination of peak probability, a probability 
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profile 300 is generated where the adjusted signaUo-noise ratio is the x-axis, 
and the probability is the y-axis. Probabflity is necessarily in the range between 
a 0% probability and a 100% probability, which is indicated as 1. Generally, the 
higher the adjusted signal-to-noise ratio, the greater the confidence that a 
5 putative peak is an actual peak. 

At some target value for the adjusted signal-to-noise/ it has been found 
that the probability is 100% that the putative peak is an actual peak and can 
confidently be used to identify the DNA composition of a biological sample. 
However, the target value of adjusted signel-to-noise ratio where the probability 

10 is assumed to be 100% is a variable parameter which is to be set according to 
application specific criteria. For example, the target signal-tonnoise ratio will be 
adjusted depending upon trial experience, sample characteristics, and the 
acceptable error tolerance in the overall system. More specifically, for situations 
requiring a conservative approach where error cannot be tolerated, the target 

16 adjusted signal-to-noise ratio can be set to, for example, 10 and higher. 

Accordingly, 100% probability will not be assigned to a peak unless the adjusted 
signal-to-nolse ratio is 10 or over. 

in other situations, a more aggressive approach may be taken as sample 
data are more pronounced or the risk of error may be reduced. In such a 

20 situation the system may be set to assume a 100% probability wfth a 6 or 
greater target signal-to-raise ratio. Of course, an intermediate sfgnai-to-nofse 
ratio target figure can be selected, such as 7, when a moderate risk of error can 
be assumed. Once the target adjusted signal-to-noise ratio is set for the method, 
then for any adjusted slgnal-to-noise ratio a probability can be determined that a 

25 putative peak is an actual peak. 

Due to the chemistry involved in performing an identification test, 
especially a mass spectrometry test of a sampie prepared by DNA amplifications, 
the allelic ratio between the signal strength of the highest peak and the signal 
strength of the second (or third and so. on) highest peak should fall within an 

30 expected ratio, if the allelic ratio falls outside of normal guidelines, the preferred 
ernbodlmerrt imposes an allelic ratio penalty to the probability. For example, HG. 
30 shows an allelic penalty 316 which has an x-axis 31 7 that is the ratio 
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between the signal strength of the second highest peak divided by signal 
strength of the highest peak. The y-axis 319 assigns a penalty between 0 and 1 
depending on the determined allelic ratio. In the preferred embodiment, it is 
assumed that allelic ratios over 30% are within the expected range and therefore 
5 no penalty is applied. Between a ratio of 10% and 30%, the penalty is linearly 
increased until at allelic ratios below 10% it is assumed the second-highest peak 
is not real. For allelic ratios between 10% and 30%, the allelic penalty chart 
315 is used to determine a penalty 319, which is multiplied by the peak 
probability determined in FIG. 29 to determine a final peak probability. Although 

10 the preferred embodiment incorporates an allelic ratio penalty to account for a 
possible chemistry error, it will be appreciated that other techniques may be 
used. Similar treatment will be applied to the other peaks. 

With the peak probability of each peak determined, the statistical 
probability for various composition components may be determined, as an 

1 5 example, in order to determine the probability of each of three possible . 
combinations of two peaks, - peak G, peak C and combinations GG, CC and 
GC. FIG. 31 shows an example where a most probable peak 325 is determined 
to have a final peak probability of 90%. Peak 325 is positioned such that it 
represents a G component in the biological sample. Accordingly, it can be 

20 maintained that there is a 90% probability that G exists in the biological sample. 
Also In the example shown in FIG. 31, the second highest probability is peak 
330 which has a peak probability of 20%. Peak 330 is at a position associated 
with a C composition. Accordingly, it can be maintained that there is a 20% 
. probability that C exists in the biological sample. 

25 With the probability of G existing (90%) and the probability of C existing 

(20%) as a starting point, the probability of combinations of G and C existing 
can be calculated. For example, FIG. 31 indicates that the probability of GG 
existing 329 is calculated as 72%. This is calculated as the probability of GG is 
equal to the probability of G existing (90%) multiplied by the probability of C not 

30 existing (100%-20%). So if the probability of G existing is 90% and the 
probability of C not existing is 80%, the probability of GG is 72%. 
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fn a similar manner, the probability of CC existing is equivalent to the 
probability of C existing (20%) multiplied by the probability of G not existing 
(100%-90%). As shown in FIG. 31, the probability of C existing is 20% while 
the probability of G not existing is 10%, so therefore the probability of CC is 
5 only 2%. Finally, the probability of GC existing is equal to the probability of G 
existing (90%) multiplied by the probability of C existing (20%), So if the 
probability of G existing is 90% and the probability of C existing is 20%, the 
probability of GC existing is 18%. In summary form, then, the probability of the 
composition of the biological sample is: 

10 probability of GG: 72%; 

probability of GC: 18%; and 
probability of CC: 2%. 
Once the probabilities of each of the possible combinations has been 
determined, FIG. 32 is used to decide whether or not sufficient confidence exists 

15 to call the genotype. FIG. 32 shows a call chart 335 which has an x-axis 337 
which is the ratio of the highest combination probability to the second highest 
combination probability. The y-axis 339 simply indicates whether the ratio is 
sufficiently high to justify calling the genotype. The value of the ratio may be 
indicated by M. The value of M is set depending upon trial data, sample 

20 composition, and the ability to accept error. For example, the value M may be 
set relatively high, such as to a value 4 so that the highest probability must be at 
least four times greater than the second highest probability before confidence is 
established to call a genotype. However, if a certain level of error may be 
- acceptable, the value of M may be set to a more aggressive value, such as to 3, 

25 so that the ratio between the highest and second highest probabilities needs to 
be only a ratio of 3 or higher. Of course, moderate value may be selected for M 
when a moderate risk can be accepted. Using the example of FIG. 31, where . 
the probability of GG was 72% and the probability of GC was 18%, the ratio 
between 72% and 18% is 4.0; therefore, whether M is set to 3, 3.5, or 4, the 

30 system would call the genotype as GG. Although the preferred embodiment 
uses a ratio between the two highest peak probabilities to determine if a 
genotype confidently can be called, it will be appreciated that other methods 
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may be substituted. It will also be appreciated that the above techniques may 
be used for calculating probabilities and choosing genotypes {or more general 
DNA patterns) consisting of combinations of more than two peaks. 

Referring now to FIG. 33, a flow chart is shown generally defining the 
5 process of statistically calling genotype described above. In FIG. 33 block 402 
shows that the height of each peak is determined and that in block 404 a noise 
profile is extrapolated for each peak. The signal is determined from the height of 
each peak in block 402 and the noise for each peak is determined using the 
noise profile in block 406. In block 410, the signal-to-noise ratio is calculated 
10 for each peak. To account for a non-Gaussian peak shape, a residual error is 
determined in block 412 and an adjusted signal-to-noise ratio is calculated in 
block 414. Block 41 6 shows that a probability profile is developed, with the 
probability of each peak existing found in block 418. An allelic penalty may be 
applied in block 420, with the allelic penalty applied to the adjusted peak 
15 probability in block 422. The probability of each combination of components is 
calculated in block 424 with the ratio between the two highest probabilities 
being determined in block 426. If the ratio of probabilities exceeds a threshold 
value, then the genotype is called in block 428. 

One skilled in the art will appreciate that processess, apparatus and 
20 systems can be practiced by other than the preferred embodiments that are 

presented in this description for purposes of illustration and not of limitation, and 
the present invention is limited only by the claims which follow. It is noted that 
equivalents for the particular embodiments discussed in this description may 
practice the invention as well. 
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WHAT IS CLAIMED IS: 

1. A method for identifying a biological sample, comprising: 

generating a data set indicative of the composition of the biological 

sample; 

5 denoising the data set to generate denoised data; 

deleting the baseline from the denoised data to generate an 
intermediate data set- 
defining putative peaks for the biological sample; 
using the putative peaks to generate a residual baseline; 
10 removing the residual baseline from the intermediate data set to 

generate a corrected data set; 

locating, responsive to removing the residual baseline, a probable 
peak in the corrected data set; and 

identifying, using the located probable peak, the biological sample. 
15 2. The method of claim 1 , wherein the data set is a spectrometry 

data set. 

3. The method of claim 1 , wherein the data set is generated by a 
mass spectrometer. 

4. The method of claim 1 , wherein denoising the data set includes 
20 generating a noise profile for the data set. 

5. The method of claim 1 , wherein denoising the data set includes 
transforming the data set using wavelet technology into a series of stages. 

6. The method of claim 5, further including generating a noise profile 
for stage 0. 

25 7. The method of claim 6, further including generating a noise profile 

for other stages. 

8. The method of claim 7, wherein the noise profile for each of the 
other stages is the noise profile for stage 0 scaled by a scaling factor. 

9. The method of claim 8, wherein the scaling factor is derived from 
30 the end portion of each of the other stages, respectively. 

10. The method of claim 5, further including applying a threshold to 
selected stages, the threshold being derived from the noise profile. 
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1 1 . The method of claim 10, wherein the threshold is scaled by a 
threshold factor before being applied to the selected stages. 

12. The method of claim 7, wherein the threshold factor is selected so 
that higher stages of data are filtered less than lower stages. 

5 13. The method of claim 5 f further including generating a sparse data 

set indicative of the denoised data. 

14. The method of claim 5, further including shifting the denoised data 
to account for variations due to a starting value for the wavelet transformation. 

15. The method of claim 1, wherein correcting the baseline further 
10 includes generating a moving average of the denoised data set. 

16. The method of claim 15, wherein the moving average is used to 
find peak sections in the denoised data set. 

17. The method of claim 16, wherein the peak sections are removed 
from the denoised data set. 

15 18. The method of claim 17, further including generating a baseline 

correction. 

19. The method of claim 1, further including compressing the 
intermediate data set, the intermediate data set having a plurality of data values 
associated with respective addresses. 
20 20. The method of claim 1 9, wherein a compressed data value is a 

real number that includes a whole portion representing the difference between 
two addresses. 

21 . The method of claim 19, wherein a compressed data value is a 
real number that includes a decimal portion representing the difference between 

25 a maximum value of all the data values and a value at a particular address. 

22. The method of claim 1, further including performing a mass shift- 
based on the position of the putative peaks. 

23. The method of claim 1, wherein generating the residual baseline 
includes deleting an area around each peak in the intermediate data. 

30 24. The method of claim 23, wherein the area deleted is derived from 

a determined width of a peak. 
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25. The method of claim 23, wherein the residual baseline is derived 
from data remaining in the intermediate data after the peaks have been removed. 

26. The method of claim 23, wherein generating the residual baseline 
includes fitting a quartic polynomial to the data remaining in the intermediate 

5 data after the peaks have been removed. 

27. The method of claim 1, wherein the probable peak is located by 
fitting a Gaussian curve to a peak area in the corrected data set. 

28. The method of claim 1 , wherein the identifying step includes using 
a generated noise profile to calculate the signal-to-noise ratio for the probable 

10 peak. 

29. The method of claim 28, wherein a residua! peak error is 
calculated by comparing the probable peak to a Gaussian curve. 

30. The method of claim 29, wherein the residual peak error is used to 
adjust the signal-to-noise ratio to generate an adjusted signal-to-noise ratio. 

15 31. The method of claim 1, wherein the identifying step includes 

deriving a peak probability for the probable peak. 

32. The method of claim 31 , wherein the peak probability is derived 
using the signal-to-noise ratio. 

33. The method of claim 31, wherein the peak probability is derived by 
20 using an allelic ratio, the allelic ratio being a comparison of two peak heights 

indicated in the corrected data. 

34. The method of claim 1 , wherein the identifying step includes 
calculating a peak probability that a probable peak in the corrected data is a peak 
indicating composition of the biological sample. 

25 35. The method of claim 34, wherein peak probability is calculated for 

each of a plurality of probable peaks in the corrected data. 

36. The method of claim 35, wherein a highest probability is compared 
to a second-highest probability to generate a calling ratio. , 

37. The method of claim 36, wherein the calling ratio is used to 
30 determine if the composition of the biological sample will be called. 

38. A system for identifying a biological sample, the system 
comprising: 



WO 02/25567 



PCTAJS01/29290 



-25- 



an instrument receiving the biological sample and generating a 
data set indicative of the composition of the biological sample; 

a computer communicating to the instrument and configured to 
receive the generated data set, the computer performing the method of: 
5 denoising the data set to generate denoised data; 

deleting the baseline from the denoised data to generate an 
intermediate data set; 

defining putative peaks for the biological sample; 
using the putative peaks to generate a residual baseline; 
10 removing the residual baseline from the intermediate data 

set to generate a corrected data set; 

locating, responsive to removing the residual baseline, a 
probable peak in the corrected data set; and 

identifying, using the located probable peak, the biological 

15 sample. 

39. The system of claim 38, wherein the computer is integral to the 
instrument. 

40. A machine readable program operating on a computing device, the 
computing device being configured to receive a data set indicating composition 

20 of a biological sample, the program implementing the steps of: 
denoising the data set to generate denoised data; 
deleting the baseline from the denoised data to generate an 
intermediate data set; 

defining putative peaks for the biological sample; 
25 using the putative peaks to generate a residual baseline; 

removing the residual baseline from the intermediate data set to 
generate a corrected data set; 

locating, responsive to removing the residual baseline, a probable 
peak in the corrected data set; and 
30 identifying, using the located probable peak, the biological sample. 

41 . A system for identifying a component of a DNA sample, 
comprising: 
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a mass spectrometer receiving the DNA sample and generating a 
data set indicative of the composition of the DNA sample; 

a computing device configured to receive the data set, the 
computing device implementing the method comprising: 
5 denoising the data set to generate denoised data: 

removing sufficiently the baseline from the denoised data to 
generate a corrected data set; 

locating a probable peak in the corrected data set; and 

identifying, using the located probable peak, a component in the 
10 composition of the DNA sample. 

42. The system of claim 41, where the method further includes using 
a statistical methodology to determine if the located probable peak is an actual 
peak. 

43. The system of claim 41 , where the method further includes 
15 determining whether the probability of the actual peak existing is sufficiently 

high to call the component of the DNA sample, and if the probability is not 
sufficiently high, then the method does not call the component. 

44. The system of claim 43, where the percentage of correctly called 
components is about 100 percent. 

20 45. A system for identifying a component in a biological sample, 

comprising: 

an instrument receiving the biological sample and generating a 
data set indicative of the component in the biological sample; 

a computing device receiving the data set and performing the 

25 steps of: 

generating corrected data by processing the data set to remove 
noise due to system and chemical reaction characteristics, the corrected data set 
having putative peak areas; 

defining the position of expected peaks using known possible peak 
30 areas from the biological sample; 

shifting the corrected data set to more closely align the putative 
peaks to the expected peaks; 
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calculating the probability that the putative peaks in the shifted 
data set are actual peaks; 

calling the composition of the biological sample responsive to the 
calculated probability. 
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